function [z]=psnr(x,y)
z=10*log10(length(x)*(32768^2)/sum((double(x)-double(y)).^2));